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Abstract 

The onset of quartetting, i.e. Q-particle condensation, in symmetric nuclear matter is studied 
with the help of an in-medium modified four nucleon equation. It is found that at very low 
density quartetting wins over pairing, because of the strong binding of the a-particles. The critical 
temperature can reach values up to around 6 MeV. Also the disappearance of a-particles with 
increasing density, i.e. the Mott transition, is investigated. In finite nuclei the Hoyle state, that is 
the O2+ of i^C, is identified as an "a-particle condensate" state. It is conjectured that such states 
also exist in heavier na-nuclei, like ^^O, ^'^Ne, etc. For instance the 6-th 0+ state of ^^O at 15.1 
MeV is identified from a theoretical analysis as being a strong candidate for an a condensate state. 
Exploratory calculations are performed for the density dependence of the a condensate fraction at 
zero temperature to address the suppression of the four-particle condensate below nuclear-matter 
density. Possible quartet condensation in other systems is discussed briefly. 

Keywords: nuclear matter, a-matter, superfluidity, Bose-Einstein condensation, strongly coupled 
systems 



1 Introduction 

One of the most amazing phenomena in quantum many-particle systems is the formation of quantum 
condensates. At present, the formation of condensates is of particular interest in strongly coupled 
fermion systems in which the crossover from Bardeen-Cooper-Schrieffer (BCS) pairing to Bose-Einstein 
condensation (BEG) may be investigated. Among very different quantum systems such as the electron- 
hole-exciton system in excited semiconductors, atoms in traps at extremely low temperatures, etc., 
nuclear matter is especially well suited for the study of correlation effects in a quantum liquid. 

Neutron matter, nuclear matter, but also finite nuclei are superfluid. However, at low density, nuclear 
matter will not cluster into pairs, i.e. deuterons but rather into a -particles which are much more stable. 
Also heavier clusters, starting with Carbon, may be of importance but are presently not considered for 
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condensation phenomena. Therefore, one may ask the question whether there exists quartetting, i.e. 
a-particle condensation, in nuclei, analogous to nuclear pairing. The only nucleus which in its ground 
state has a pronounced a -cluster structure is ^Be. In section 4 we will show a figure of ^Be in the 
laboratory frame and in the intrinsic deformed frame. We will see that ^Be is formed out of two a- 
particles roughly 4 fm apart, only mildly interpenetrating one another. Actually ^Be is slightly unstable 
and the two a's only hold together via the Coulomb barrier. Because of the large distance of the two 
a-particles, the 0^ ground state of ^Be has, in the laboratory frame, a spherical density distribution 
whose average is very low: about 1/3 of ordinary saturation density po- ^Be is, therefore, a very large 
object with an rms radius of about 3.7 fm to be compared with the nuclear systematics oi R = tqA^^^ = 
2.44 fm. Definitely ^Be is a rather unusual and, in its kind, unique nucleus. One may ask the question 
what happens when one brings a third a-particle alongside of ^Be. We know the answer: the 3-a system 
collapses to the ground state of ^^C which is much denser than ^Be and can not accommodate with its 
small radius of 2.4 fm three more or less free a-particles barely touching one another. One nevertheless 
may ask the question whether the dilute three a configuration ^Be-a, or rather a-a-a, may not form 
an isomeric or excited state of ^^C. That such a state indeed exists will be one of the main subjects 
of our considerations. Once one accepts the idea of the existence of an a-gas state in ^^C, there is no 
reason why equivalent states at low density should not also exist in heavier na- nuclei, like ^^O, ^°Ne, 
etc. The possible existence of a loosely bound Aa state in ^^O will be another topic of our presentation. 
In a mean field picture, i.e. all a's being ideal bosons ( in this context remember that the first excited 
state of an a-particle is at ~ 20 MeV, by factors higher than in all other nuclei), all a's will occupy 
the lowest OS-state, i.e. they will condense. This forms, of course, not a macroscopic condensate but it 
can be understood in the same sense as we know that nuclei are superfiuid because of the presence of a 
finite number of Cooper pairs. On the other hand, for example during the cooling process of compact 
stars pp, where one predicts the presence of a-particles [2], a real macroscopic phase of condensed a's 
may be formed. In the present contribution we will mainly concentrate on nuclear systems but we also 
can think about the possibility of quartetting in other Fermi-systems. One should, however, keep in 
mind that a pre-requisite for its existence is, as in nuclear physics, that there are four different types of 
fermions. For example to form quartets with cold atoms one could try to trap fermions in four different 
magnetic substates, a task which eventually seems possible [5]. Also theoretical works in this direction 
have appeared in the mean while [H [S] . The fact that the a-particle condensates in nuclei do not form 
the ground state, may give raise to questions. In this respect, one should note that Bose condensates 
of cold atoms in traps also are not in their ground state which is a solid. It is a question of time scales, 
a-particle condensate states in nuclei usually live four orders of magnitude longer than typical nuclear 
times. 

In the next section we will investigate how the binding energy of various nuclear clusters change 
with density. In section 3 we study the critical temperature of a-particle condensation in infinite 
matter via an in-medium four-nucleon equation (Thouless criterion) and in section 4 we treat a-particle 
condensation in finite nuclei. In section 5 we give our results for ^^C and in section 6 the ones for 
^^O. In section 7, we briefiy discuss the question of the occupation numbers of the a-particles in the 
various nuclear states with emphasis on the condensate states. In section 8 we present a simplified 
calculation of the condensate fraction in a-matter. Finally in section 9 we conclude with outlook and 
further discussions. 

2 Nuclear clusters in the medium 

With increasing density of nuclear matter, medium modifications of single-particle states as well as of 
few-nucleon states become of importance. The self-energy of an A-particle cluster can in principle be 
deduced from contributions describing the single-particle self-energies as well as medium modifications of 
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the interaction and the vertices. A guiding principle in incorporating medium effects is the construction 
of consistent ("conserving") approximations, which treat medium corrections in the self-energy and in 
the interaction vertex at the same level of accuracy. This can be achieved in a systematic way using 
the Green functions formalism [6]. At the mean-field level, we have only the Hartree-Fock self-energy 
pHF _ V{12, 12)ex/(2) together with the Pauli blocking factors, which modify the interaction from 
V{12, 1'2') to 1/(12, 1'2')[1 - /(I) - /(2)], with /(I) = [1 + exp{E^^{l) - ii)/T]-\ In the case of the 
two-nucleon system {A = 2), the resulting effective wave equation which includes those corrections reads 

[^^^(1) + ^^^(2) - E2,n,p] V^2,n,p(12) + ^[1 - /(I) - /(2)] V^(12, l'2')V^2,n,p(l'2') = 0. (1) 

1'2' 

This effective wave equation describes bound states as well as scattering states. The onset of pair 
condensation is achieved when the binding energy Ed,p=o coincides with 2/i. 

Similar equations have been derived from the Green function approach for the case A = 3 and 
A = 4, describing triton/helion (^He) nuclei as well as a-particles in nuclear matter. The effective wave 
equation contains in mean field approximation the Hartree-Fock self-energy shift of the single-particle 
energies as well as the Pauli blocking of the interaction. We give the effective wave equation for A = 4, 



+ B«f (2) + B«"'(3) + B«f (4) - B4.„.p] '*4,„.p(1234) 



i<j 1'2'3'4' 



fijWitj,t'j') n 4,fc'^4,n,p(l'2'3'4') = 0. 



(2) 



A similar equation is obtained for A = 3. 

The effective wave equation has been solved using separable potentials for A = 2 hj integration. 
For A = 3,4 we can use a Faddeev approach The shifts of binding energy can also be calculated 
approximately via perturbation theory. In Fig. [1] we show the shift of the binding energy of the light 
clusters {d,t/h and a) in symmetric nuclear matter as a function of density for temperature T = 10 
MeV. 




Figure 1: Shift of binding energy of the light 
clusters {d - dash dotted, t/h - dotted, and a 
- dashed: perturbation theory, full line : non- 
perturbative Faddeev- Yakubovski equation) in 
symmetric nuclear matter function of den- 
sity for given temperature T = 10 MeV [7j. 



nuclear density p[fm ] 



It is found that the cluster binding energy decreases with increasing density. Finally, at the Mott 
density p^°n^p{T) the bound state is dissolved. The clusters are not present at higher densities, merging 
into the nucleonic medium. For a given cluster type characterized by A, n, we can also introduce the 
Mott momentum P^'^^{p,T) in terms of the ambient temperature T and nucleon density p, such that 
the bound states exist only for P > P^'^^{p, T). We do not present an example here, but it is intuitively 
clear that a cluster with high c.o.m. momentum with respect to the medium is less affected by the 
Pauli principle than a cluster at rest. 
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chemical potential n* [MeV] density n, [fm ^] 



Figure 2: Transition temperature to quar- 
tetting/ pairing function of chemical 
potential in symmetric nuclear matter. 



Figure 3: Transition temperature to quartet- 
ting/pairing as a function of nucleon density in 
symmetric nuclear matter. 



3 Four-particle condensates and quartetting in nuclear mat- 
ter 

In general, it is necessary to take into account of all bosonic clusters to gain a complete picture of the 
onset of superfluidity. As is well known, the deuteron is weakly bound as compared to other nuclei. 
Higher A-clusters can arise that are more stable. In this section, we will consider the formation of a- 
particles, which are of special importance because of their large binding energy per nucleon (~ 7 MeV). 
We will not include tritons or helions, which are fermions and not so tightly bound. Moreover, we 
will not consider nuclei in the iron region, which have even larger binding energy per nucleon than the 
a-particle and thus constitute, in principle, the dominant component at low temperatures and densities. 
However, the latter are complex structures of many particles and are strongly affected by the medium 
as the density increases for given temperature, so that they are assumed not to be of relevance in the 
density region considered here. 

The in-medium wave equation for the four-nucleon problem has been solved using the Faddeev- 
Yakubovski technique, with the inclusion of Pauli blocking. The binding energy of an a-like cluster 
with zero c.o.m. momentum vanishes at around po/10, where po — 0.16 nucleons/fm^ denotes the 
saturation density of isospin-symmetric nuclear matter, see Fig. [1] Thus, the four-body bound states 
make no significant contribution to the composition of the system above this density. Given the medium- 
modified bound-state energy -E'4,p, the bound-state contribution to the EOS is 

P 

We will not include the contribution of the excited states or that of scattering states. Because of the 
large specific binding energy of the a particle, low-density nuclear matter is predominantly composed 
of a particles. This observation underlies the concept of a matter and its relevance to diverse nuclear 
phenomena. 

As exemplified by Eq. ([2]), the effect of the medium on the properties of an a particle in mean- 
field approximation (i.e., for an uncorrelated medium) is produced by the Hartree-Fock self-energy 
shift and Pauli blocking. The shift of the a-like bound state has been calculated using perturbation 
theory [8] as well as by solution of the Faddeev-Yakubovski equation [Tj. It is found that the bound 
states of clusters d, t, and h with A < 4 are already dissolved at a Mott density ~ Po/10, 

see Fig. [H Since Bose condensation only is of relevance for d and a, and the fraction of d, t and 
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h becomes low compared with that of a with increasing density, we can neglect the contribution of 
them to an equation of state. Consequently, if we further neglect the contribution of the four-particle 
scattering phase shifts in the different channels, we can now construct an equation of state p(T, /i) = 
p^^'^'^iT, fi) + p^°'^^'^''^(T, fi) + p^°'^^'^'°'(T, fi) such that a-particles determine the behavior of symmetric 
nuclear matter at densities below p^°** and temperatures below the binding energy per nucleon of the 
a-particle. The formation of deuteron clusters alone gives an incorrect description because the deuteron 
binding energy is small, and the abundance of d-clusters is small compared with that of a-clusters. In 
the low density region of the phase diagram, a-matter emerges as an adequate model for describing the 
nuclear-matter equation of state. 

With increasing density, the medium modifications - especially Pauli blocking - will lead to a 
deviation of the critical temperature Tc(p) from that of an ideal Bose gas of a-particles (the analogous 
situation holds for deuteron clusters, i.e., in the isospin-singlet channel). 

Symmetric nuclear matter is characterized by the equality of the proton and neutron chemical 
potentials, i.e., Pp = Pn = Then an extended Thouless condition based on the relation for the 
four-body T-matrix (in principle equivalent to Eq. (|2]) at eigenvalue Ap) 



T4(1234,r'2"3"4",4/i) 



^ y(12,r20[l-/(l)-/(2)] 

+cycl. iT4(l'2'3'4', l"2"3"4",4p) 



5(3,305(4,4') 



(4) 



serves to determine the onset of Bose condensation of a-like clusters, noting that the existence of a 
solution of this relation signals a divergence of the four-particle correlation function. An approximate 
solution has been obtained by a variational approach, in which the wave function is taken as Gaussian 
incorporating the correct solution for the two-particle problem [9]. 

The results are presented in Figs. [2] and [31 An important consequence of those is that at the lowest 
temperatures, Bose-Einstein condensation occurs for a particles rather than for deuterons. As the 
density increases within the low-temperature regime, the chemical potential p first reaches —7 MeV, 
where the a's Bose-condense. By contrast, Bose condensation of deuterons would not occur until p 
rises to —1.1 MeV. 

The "quartetting" transition temperature sharply drops as the rising density approaches the critical 
Mott value at which the four-body bound states disappear. At that point, pair formation in the isospin- 
singlet deuteron-like channel comes into play, and a deuteron condensate will exist below the critical 
temperature for BCS pairing up to densities above the nuclear-matter saturation density po, as described 
in the previous Section. The critical density at which the a condensate disappears is estimated to be 
Po/3. Therefore, a-particle condensation primarily only exists in the Bose- Einstein-Condensed (BEC) 
phase and there does not seem to exist a phase where the quartets acquire a large extension as Cooper 
pairs do in the weak coupling regime. However, the variational approach of Ref. [H] on which this 
conclusion is based represents only a first attempt at the description of the transition from quartetting 
to pairing. The detailed nature of this fascinating transition remains to be clarified. 

Many different questions arise in relation to the possible physical occurrence and experimental 
manifestations of quartetting: Can we observe the hypothetical "a condensate" in nature? What about 
thermodynamic stability? What happens with quartetting in asymmetric nuclear matter? Are more 
complex quantum condensates possible? What is their relevance for finite nuclei? As discussed below, 
the special type of microscopic quantum correlation associated with quartetting may be important in 
nuclei, its role in these finite inhomogeneous systems being similar to that of pairing. 
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4 Description of Alpha-Particle Condensate States in Self- 
Conjugate 4n Nuclei 



Figure 4: Contours of constant density (taken from 
Ref. [in]), plotted in cylindrical coordinates, for 
®Be(0~''). The left side (a) is in the "laboratory" frame 
while the right side (b) is in the intrinsic frame. 



"■-4 -2 2 4 -2 2 4 

t = (i2+y2)i« (fm) 

Let us discuss the possibility of quartetting in nuclei. The only nucleus having a pronounced a- 
cluster structure in its ground state is ®Be. In Fig. |4](a), we show the result of an exact calculation of 
the density distribution of ®Be in the laboratory frame. In Fig. Hl^b) we show, for comparison, the result 
of the same calculation in the intrinsic, deformed frame. We see a pronounced two a-cluster structure 
where the two a's are ~ 4 fm apart, giving rise to a very low average density p ~ Po/3 as seen in 
Fig. m^a). As already discussed in the introduction, *Be is a rather unusual and unique nucleus. One 
may be intrigued by the question, already raised earlier, whether loosely bound a-particle configurations 
may not also exist in heavier na-nuclei, at least in excited states, naturally close to the na disintegration 
threshold. Since a-particles are rather inert bosons ( first excited state at ~ 20 MeV), these a-particles 
then would all condense in the lowest S-wavefunction, very much in the same way as do bosonic atoms 
in magneto-optical traps [TTj. This question and exploring related issues of quartetting in finite nuclei 
will consume most of the rest of the present study. 

In fact, we will be able to offer strong arguments that the 0^ state of ^^C at 7.654 MeV is a state 
of a-condensate nature. Later, also indications for the existence of an analogous state in ^^O will be 
discussed. 

First, it should be understood that the 0^ state in ^^C is in fact hadronically unstable (as *Be), 
being situated about 300 keV above the three a-break up threshold. This state is stabilized only by the 
Coulomb barrier. It has a width of 8.7 eV and a corresponding lifetime of 7.6 x 10~^^ s. As well known, 
this state is of paramount astrophysical (and biological!) importance due to its role in the creation 
of ^^C in stellar nucleosynthesis. Its existence was predicted in 1953 by the astrophysicist Fred Hoyle 
[12]. His prediction was confirmed experimentally a few years later by Willy Fowler and coworkers at 
Caltech [13]. It is also well known that this Hoyle state, as it is now called, is a notoriously difficult 
state for any nuclear theory to explain. For example, the most modern no-core shell-model calculations 
predict the 0^ state in ^^C to lie at around 17 MeV above the ground state - more than twice the 
actual value [H]. This fact alone tells us that the Hoyle state must have a very unusual structure. It is 
easy to understand that, should it indeed have the proposed loosely bound three a-particle structure, 
a shell-model type of calculation would have great difficulties in reproducing its properties. 

An important development bearing on this issue took place some thirty years ago. Two Japanese 
physicists, M. Kamimura [T3] and K. Uegaki [TB], along with their collaborators, almost simultaneously 
reproduced the Hoyle state from a microscopic theory. They employed a twelve-nucleon wave function 
together with a Hamiltonian containing an effective nucleon-nucleon interaction. At that time, their 
work did not attract the attention it deserved; the true importance of their achievement has been 
appreciated only recently. The two groups started from practically the same ansatz for the ^^C wave 
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function, which has the following three a-cluster structure: 

(ri...ri2r'C)=^[x(s,t)0i0203] . (5) 

In this expression, the operator A imposes antisymmetry in the nucleonic degrees of freedom and (pi, 
with i = 1,2,3 for the three a's, is an intrinsic a-particle wave function of prescribed Gaussian form, 

(/.(ri,r2,r3,r4) =exp{- [(n - ra)^ + (n - rg)^ + ■ ■ ■ ] /Sfe^} , (6) 

where the size parameter b is adjusted to fit the rms of the free a-particle, and is a yet-to-be 

determined three-body wave function for the c.o.m. motion of the three «'s, their corresponding Jacobi 
coordinates being denoted by s and t. The unknown function x was determined via calculations based 
on the Generator Coordinate Method [16] (GCM) and the Resonating Group Method [15] (RGM) 
calculations using the Volkov I and Volkov II nucleon-nucleon forces, which fit a-a phase shifts. The 
precise solution of this complicated three body problem, carried out three decades ago, was truly a 
pioneering achievement, with results fulfilling expectations. The position of the Hoyle state, as well as 
other properties including the inelastic form factor and transition probability, successfully reproduced 
the experimental data. Other states of ^^C below and around the energy of the Hoyle state were also 
successfully described. Moreover, it was already recognized that the three a's in the Hoyle state form 
sort of a gas-like state. In fact, this feature had previously been noted by H. Horiuchi [T7] prior to 
the appearance of Refs. [151 [IS], based on results from the orthogonality condition model (OCM) |18j . 
All three Japanese research groups concluded from their studies that the linear-chain state of three 
a-particles, postulated by Morinaga many years earlier [19] as an interpretation of the Hoyle state, had 
to be rejected. 

Although the evidence for interpreting the Hoyle state in terms of an a gas was stressed in the cited 
papers from the late 1970's, two important aspects of the situation were missed at that time. First, 
because the three a's move in identical S-wave orbits, one is dealing with an a-condensate state, albeit 
not in the macroscopic sense, and that this may be a quite general phenomenon, also in heavier self- 
conjugate nuclei. The second important point is that the complicated three-body wave function x(s, t) 
for the c.o.m. motion of the three a's can be replaced by a structurally and conceptually very simple 
microscopic three-a wave function of the condensate type, which has practically 100 percent overlap 
with the previously constructed ones [20] [21] (see also Ref. [22]). We now describe this condensate 
wave function. 

We start by examining the BCS wave function of ordinary fermion pairing, obtained by projecting 
the familiar BCS ground-state ansatz onto an A^-particle subspace of Fock space. In the position 
representation, this wave function is 

(ri ■ ■ ■ r^,|BCS) = A [0(ri, r^mr-,, r,) ■ ■ ■ <P{r^_^r^)] , (7) 

where is the Cooper-pair wave function (including spin and isospin), which is to be determined 

variationally through the familiar BCS equations. The condensate character of the BCS ansatz is borne 
out by the fact that within the antisymmetrizer A, one has a product of A^/2 times the same pair wave 
function 0, with one such function for each distinct pair in the reference partition of {1, 2, . . . , iV — 1, N}. 
Formally, it is now a simple matter to generalize ([7]) to quartet or a-particle condensation. We write 

(ri, . . . , VNl^na) = -A. [0a(^^l, ^^2, '''S, 'r^(\)^{r^, . . . , rg) ■ ■ ■ 0a('^Af-3, • • • , 'Tn)] , (8) 

where (pa is the wave function common to all condensed a-particles. Of course, finding the variational 
solution for this function is, in general, extraordinarily more complicated than finding the Cooper pair- 
wave function of Eq. ([7]). Even so, in the present case that the a-particle is the four-body cluster 
involved, and for applications to relatively light nuclei, the complexity of the problem can be reduced 
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dramatically. This possibility stems from the fact, already known to the authors of Refs. [151 [16], that, 
due to the BEC-character of a-particle condensation (see above), an excellent variational ansatz for 
the intrinsic wave function of the a-particle is provided [as in Eq. ([H])], by a Gaussian form with only 
the size parameter b to be determined. In addition - and here resides the essential point of our wave 
function - even the c.o.m. motion of the system of a-particles can be described very well by a Gaussian 
wave function with, this time, a size parameter B ^ b to account for the motion over the nuclear space. 
We therefore write ^ ^ 

Mri, r2, rs, r^) = e~'^^' I ^ 0(ri - ra, ri - rg, ■ ■ ■ ) , (9) 

where -R = (^i + ^2 + + r4)/4 is the c.o.m. coordinate of one a-particle and 0(ri — t'2, ...) is the 
same intrinsic a-particle wave function of Gaussian form as already used in Refs. [T5l [T6] and given 
explicitly in Eq. ([6]). Naturally, in Eq. ([8]) the center of mass of the three a's, i.e., of the whole 
nucleus, should be eliminated; this is easily achieved by replacing Rhj R — X^m in each of the a wave 
functions in Eq. ([H]). The a-particle condensate wave function specified by Eqs. ((HI) and Qj, proposed 
in Ref. [20] and henceforth called the THSR wave function, now depends on only two parameters, B 
and b. The wave function ([8]) with ([9]) is pictorially represented in Fig. [51 The expectation value of an 
assumed microscopic Hamiltonian H, 



n{B,b) 



{<^naiB,b)\H\<!>naiB,b)) 



(10) 



can be evaluated, and the corresponding two-dimensional energy surface can be quantized using the 
two parameters B and b as Hill- Wheeler coordinates. 




B 



Figure 5: Pictorial representation of the THSR wave function for n = 3 (^^C). The three a-particles are 
trapped in the OS-state of a wide harmonic oscillator (B) and the four nucleons of each a are confined 
in the Os-state of a narrow one (6). All nucleons are antisymmetrised. 



Before presenting the results, let us discuss the THSR wave function in somewhat more detail. This 
innocuous- looking variational ansatz, namely Eq. ([8]) together with Eq. (Q, is actually more subtle than 
it might at first appear. One should realize that two limits are incorporated exactly. One is obtained 
by choosing B = b, for which Eq. ([8]) reduces to a standard Slater determinant with harmonic-oscillator 
single-nucleon wave functions, leaving the oscillator length b as the single adjustable parameter. This 
holds because the right-hand-side of expression with B = b, becomes a product of four identical 
Gaussians, and the antisymmetrization creates all the necessary P, D, etc. harmonic oscillator wave 
functions automatically [20]. On the other hand, when B ^ b, the density of a-particles is very low, 
and in the limit B ^ oo, the average distance between a's is so large that the antisymmetrisation 
between them can be neglected, i.e., the operator A in front of Eq. ([HD becomes irrelevant and can be 
removed. In this limiting case, our wave function then describes an ideal gas of independent, condensed 
a-particles - it is a pure product state of a's! An elucidating study on this aspect is given in Ref. [23] . 
Evidently, however, in realistic cases the antisymmetrizer A cannot be neglected, and evaluation of the 
expectation value (fTOl) becomes a nontrivial analytical task. The Hamiltonian in Eq. (ITOl) was taken to 
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be the one used in Ref. [2l] , which features an effective nucleon-nucleon force of the Gogny type, with 
parameters fitted to a-a scattering phase shifts as available about fifteen years ago. This force also leads 
to very reasonable properties of ordinary nuclear matter. Our theory is therefore free of any adjustable 
parameters. The energy landscapes T-C{B,b) for various na nuclei are shown in Fig. [6] [251 ISSj- We 
see that they all have qualitatively the same structure. From the minimum point on, with increasing 
5-parameter, developes a valley with constant b which takes approximately the value of the free space 
a-particle. The valley then goes over a saddle point which indicates the disintegration of the nucleus 
into n a-particles. It is interesting to see that the minimum of the energy surfaces does not correspond 
to the Slater determinant case with b = B but rather a quite substantial gain in energy due to four 
body correlations can be observed even for the ground state in these light nuclei. 




Figure 6: Contour map of the energy surface Ti.{B, b) for (a) ®Be, (b) ^^C and (c) ^^O. The variables B 
and i?o are connected by the relation, B'^ = b'^ + 2Rq. Numbers attached to the contour lines are the 
binding energies in units of MeV. 



It is evident that for large numbers of a-particles the explicit antisymmetrisation demanded in our 
wave function of Eq. ([8]), will encounter great difficulties. On the other hand in the case of many 
quartets or a macroscopic number of it in the condensate, in particular in infinite nuclear matter, one 
may transform our number conserving condensate wave function into a coherent state and proceed in a 
similar way as in the BCS case of standard pairing. Roughly this can go along the following lines. The 
first step is to define an a-particle coherent state 
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|«) ~ exp (^^ 52344444) |vac). (11) 

■ 1234 

Next comes to mix pair creators and pair destructors: 

Qi = l E{^"244 - ^12ClC2}. (12) 
12 

If we replaced the pair operators by bosons, this would correspond to a Hartree-Fock-Bogoliubov 
approach for bosons, with, in analogy to the fermion case, a bosonic gap equation, etc. However, here 
we want to keep the fermionic structure of the pairs. The amplitudes X, Y shall obey the following 
orthonormality, respectively completeness relations XX'^ — YY'^ = 1. We also demand: Qa\<y) = 0, and 
the corresponding order parameter is a quartet expectation value: (c|444) ~ XY. 

A self-consistent set of quartet equations can then be constructed: 

^ ^'^^-^^ e(^.), (13) 



-A ) \ Y ) \ Y 



where 

A4 = ([ctct, [H, ctct]]) ~ v„M4clcl), (14) 

and f . .. stands for the matrix element of the interaction and the bosonic "gap" is then very schematically 
given by 

A4 ~ ^t;....XF. (15) 

Similar equations from an analogous procedure can be obtained mixing a fermion destructor with 
three creators: gj^ = X] ^123^144 ~ Xl^i""^!- The final equation for quartetting is a very intuitive 
extension of the equation for the BCS pairing order parameter (cc) 



(ei + €2 + 63 + eAK'' 



1234 



Til — 722) A°234 + perms. 



A^K^,,,, (16) 



with ej the self-consistent single particle energies and A"234 = \ 'Vi22,'i'K^ii>zi where -ft'1234 is the 
order parameter specified below, and W1234 is the antisymmetrised matrix element of the interaction. 
The occupation numbers rii = {c\ci) are obtained from the single particle Green's function Gi = 
{u-ei-M'^y^ with the mass operator M'^ = ^ ^234 1^234^^234 (^+e2+e3+e4)"^ with A^^^ = N^s^D^u 
and Ai'234 = (1 — n2)(l — 7^3) (1 — n^) + n2n3n4 (all energies counted from /i). In the zero density limit 
Eq. (fT6l) goes over into the exact free space a-particle Schrodinger equation. 

In infinite matter, the order parameter has zero total momentum and we write 

i^l234 = (4444) ^ <^(fcl + ^2 + ^3 + k^){('k/k/k/k}sT ^^^^ 

where S, T stand for the spin-isospin wave function. Thus, the order parameter depends on three 
momenta n = ki — k2; k' = k^ — k^, P = ki + k2 = —{k^ + ^4). In general this makes 9 variables! 
Nonlinear equations for the order parameter have to be solved what is a very demanding numerical 
task. 

However, strong simplifications may be possible! For Bose-Einstein condensation of fermion-clusters 
heavier than a pair, one may proceed to a mean-field description of the cluster with projection on good 
total momentum K. The clusters condense in the K = state. Even for the a-particle a mean field 
description is quite a good approximation, under the condition to use effective forces, such as Skyrme 
or Gogny, and to project, as mentioned, on good linear momentum (see Fig. [7]). 

~ 6{ki + k2 + k3 + ki)ifosiki)iposik2)(posik3)(fosik4). (18) 
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Figure 7: Spin-isospin saturated (05*)^ mean field configuration of the ct-particle. 



With 



(a) 



(YX )l234, 



(19) 



1234 



and 



XX^ - YY^ = 1 XX^ = (1 - $(")^) 



-1 



(20) 



the self-consistent quartet equation can be expressed entirely in terms of and then via the product 
ansatz everything via a single OS* wavefunction. The nonlinear equation for (pos{f) should be solvable! 
This for any number of a-particles! The last part of this section is onging work with T. Sogo [28] . 




As we discussed already, the variational wave function constructed from the Hill- Wheeler equation 
based on Eqs. ([8]), Qj, and ffTOj) has practically 100 percent overlap with the RGM and GCM wave 
functions constructed in Refs. [TS] and [TB], once the same Volkov force is used [2T]. It is, therefore, not 
astonishing that our results are very similar to theirs. Nevertheless, let us again discuss the situation 
in some detail. For ^^C we obtain two eigenvalues in the 0"'"-channel: the ground state and the Hoyle 
state. Theoretical values for positions, rms values, and transition probabilities are given in Tableland 
compared to the data. Inspecting the rms radii, we see that the Hoyle state has a volume 3 to 4 larger 
than that of the ground state of ^^C. This is the primary aspect of the dilute-gas state we highlighted 
above. We also can make a deformed calculation in allowing the width parameter B to have different 
values in the different directions. Projecting on good angular momentum then yields the position of 
the second 2"'"-state in ^^C which is in good agreement with the experimental value [29l [30] . Also its 
width can be evaluated and one obtains a quite reasonable estimate. Detailed investigation of the wave 
function of the 2^-state shows that it can essentially be described in lifting out of the condensate state 
with the three a's in the OS'-orbit, one a-particle in the next Oi^-orbit. It is tempting to imagine that 
the Ojj'-state which, experimentally, is almost degenerate with the 2^-state, is obtained by lifting one 
a-particle into the IS'-orbit. Preliminary theoretical studies [31] indicate that this scenario might indeed 
apply. However, the width of the 0^ state is very broad (~ 3 MeV), rendering a theoretical treatment 
rather delicate. Further investigations are necessary to validate or reject this picture which is shown 
graphically in Fig. O At any rate, it would be quite satisfying, if the triplet of states, (0^, 2^, O3 ) could 
all be explained from the a-particle perspective, since those three states are precisely the ones which 
cannot be explained within a (no core) shell model approach |14j . 

Constructing a pure-state a-particle density matrix p{R, R') from our wave function, integrating out 
of the total density matrix all intrinsic a-particle coordinates, and diagonalizing this reduced density 
matrix, we find that the corresponding OS* a-particle orbit is occupied to 70 percent by the three a- 
particles [231 [22] • This is a huge percentage, giving vivid support to the view that the Hoyle state is an 
almost ideal a-particle condensate. We will dwell on this point in more detail in sections 7 and 8. 

Let us now discuss what to our mind is the most convincing evidence that our description of the 
Hoyle state is the correct one. Like the authors of Ref. [15], we reproduce very accurately the inelastic 
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Figure 8: (Color online) Theoretical inter- 
pretations of the 0^, 2^ and O3 states, which 
are missing in the results of no-core shell- 
model calculation (see Fig. [9]), are schemati- 
cally shown. 



Figure 9: Experimental and no core shell model 
positive-parity excitation spectra of ^^C. The 
effective interaction was derived from the CD- 
Bonn NN potential in a HO basis with Ml = 15 
MeV. Figure taken from 



form factor 0^ 0^ of ^^C, as shown in Fig. [101 As such, the agreement with experiment is already 
quite impressive. Additionally, however, the following study was made, results from which are presented 
in Fig. [m We artificially varied the extension of the Hoyle state and examined the influence on the 
form factor. It was found that the overall shape of the form factor shows little variation, for example 
in the position of the minimum. On the other hand, we found a strong dependence of the absolute 
magnitude of the form factor; Fig. [TT] illustrates this behavior with a plot showing the variation of the 
height of the first maximum of the inelastic form factor as a function of the percentage change of the 
rms radius of the Hoyle state [33]. It can be seen that a 20 percent increase of the rms radius produces 
a remarkable decrease of the maximum by a factor of two! This strong sensitivity of the magnitude 
of the form factor to the size of the Hoyle state enhances our firm belief that the agreement with the 
actual measurement is tantamount to a proof that the calculated wide extension of the Hoyle state 
corresponds to reality. 

Summarizing our inquiry into the possible role of a clustering in ^^C, we have accumulated enough 







condensate w.f. 
(Hill-Wheeler) 


RGM [15] 


Exp. 


^(MeV) 


0^ 
0+ 


-89.52 
-81.79 


-89.4 
-81.7 


-92.2 
-84.6 


i?r.m.s.(fm) 


0^ 
0+ 


2.40 
3.83 


2.40 
3.47 


2.44 


M(0+ ^0+)(fm2) 




6.45 


6.7 


5.4 



Table 1: Comparison of the binding energies, rms radii (-Rr.m.s.), and monopole matrix elements 
(M(0^ — > ^i)) for ^^C given by solving Hill-Wheeler equation based on Eq. ([8]) and by Ref. |15j . 
The effective two-nucleon force Volkov No. 2 was adopted in the two cases for which the 3a threshold 
energy is calculated to be —82.04 MeV. 
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Figure 10: Experimental values of inelastic 
form factor in ^^C to the Hoyle state are 
compared with our values and those given by 
Kamimura et al. in Ref. [13] (RGM). In our 
result, the Hoyle-state wave function is ob- 
tained by solving the Hill- Wheeler equation 
based on Eq. ([8]). The experimental values 
are taken from Ref. 



Figure 11: The ratio of the value of the 
maximum height, theory versus experi- 
ment, of the inelastic form factor, i.e. 
max|F(g)p/max|F(g)|g^p, is plotted as a 
function of 6 = (-Rr.m.s. — Ro)/Ro- Here 
-Rr.m.s. and Ro are the rms radii corre- 
sponding, respectively, to the wave func- 
tion of Eq. ([H]) and that obtained by solv- 
ing the Hill- Wheeler equation based on 
Eq. iD. 



facts to be convinced that the Hoyle state is, indeed, what one may call an a-particle condensate state. 
At the same time, we acknowledge that referring to only three particles as a "condensate" constitutes 
a certain abuse of the word. However, in this regard it should be remembered that also in the case 
of nuclear Cooper pairing, only a few pairs are sufficient to obtain clear signatures of superfluidity in 
nuclei! 

Let us now go one step further and investigate the four a-particle case. 

6 Alpha-particle condensation in ^^O 

The establishment of this condensate aspect of the Hoyle state naturally leads us to the speculation 
about 4a-particle condensation in ^^O, which is the focus in this section. The situtaion in ^^O is, as 
compared to ^^C quite a bit more complicated, even in the O"*" channel alone. This stems from the fact 
that while "knocking loose" one a-particle in ^^C, necessarily the other two are also almost free (^Be) 
and all three a's form the gas state. However, exciting one a-particle out of the ground state in ^^O, 
may leave the remaining ^^C core in the ground state or in various excited states of the shell model 
type. Therefore, in ^^O we need to "knock loose" at least two a's to obtain the a gas state. 

The O"*" spectrum of ^^O has, in the past, very well been reproduced up to about 13 MeV excitation 
energy, including the ground state, with a semi-microscopic cluster model, i.e. the a+^^C OCM (Or- 
thogonality Condition Model) [35]. In particular, this model calculation, as well as that of an a+^^C 
Generator-Coordinate-Method one [36], demonstrates that the O2 state at 6.05 MeV and the O3 state 
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at 12.05 MeV have a+^'^C structures [37] where the a-particle orbits around the ^^C(0||^)-core in an 
S'-wave and around the ^^C(2^)-core in a -D-wave, respectively. Consistent results were later obtained 
by the 4a OCM calculation within the harmonic oscillator basis [SB]- However, the model space adopted 
in Refs. [351 EH [38] is not sufficient to account simultaneously for the a+^^C and the 4a gas-like config- 
urations. On the other hand, the 4a-particle condensate state was first investigated in Ref. [20] and its 
existence was predicted around the 4a threshold with the a-particle condensate wave function. While 
this so-called THSR wave function can well describe the dilute a cluster states as well as shell model 
like ground states, other structures such as a+^^C clustering are smeared out and only incorporated in 
an average way. Since there exists no calculation, so far, which reproduces both the 4a gas and a+^'^C 
cluster structures simultaneously, it is crucial to perform an extended calculation for the simultaneous 
reproduction of both kinds of structures, which will give a decisive benchmark for the existence of the 
4a-particle condensate state from a theoretical point of view. 

Therefore the objective in Ref. [3H] was to explore the 4a condensate state by solving a full OCM 
four-body equation of motion without any assumption with respect to the structure of the 4a system. 
Here we take the 4a OCM with Gaussian basis functions, the model space of which is large enough to 
cover the 4a gas, the a-|-^^C cluster, as well as the shell-model configurations. The OCM is extensively 
described in Ref. [18]. Many successful applications of OCM are reported in Ref. [ID]. The 4a OCM 
Hamiltonian is given as follows: 



4 4 4 

n = Y^T,-T,^ + Y,[vS\^,J)+V,^^\^,J)+vi^\^,J)] + J] V3o.{^,J,k) +V,^{1,2,3,4), (21) 

i i<j i<j<k 

where Tj, V2a\i,j), y2a\hj)^ V3a{i,j,k) and V4a(l,2,3,4) stand for the operators of kinetic energy 
for the i-th a particle, two-body. Coulomb, three-body and four-body forces between a particles, re- 
spectively. The center-of- mass kinetic energy is subtracted from the Hamiltonian. V2^\i,j) is the 
Pauli exclusion operator [41], by which the Pauli forbidden states between two a-particles in OS", OD 
and IS* states are eliminated, so that the ground state with the shell-model-like configuration can be 
described correctly. The effective a-a interaction V2a^ is constructed by the folding procedure from two 
kinds of effective two-nucleon forces. One is the Modified Hasegawa-Nagata (MHN) force [12] and the 
other is the Schmidt- Wildermuth (SW) force [13], see Refs. [23] and [31] for applications, respectively. 
We should note that the folded a-a potentials reproduce the a-a scattering phase shifts and energies 
of the ^Be ground state and of the Hoyle state. The three-body force Vsq, is as in Refs. [23] and [3Tj 
where it was phenomenologically introduced, so as to fit the ground state energy of ^^C. In addition, 
the phenomeno logical four-body force V^a which is taken to be a Gaussian is adjusted to the ground 
state energy of ^^O, where the range is simply chosen to be the same as that of the three-body force. 
The origin of the three- body and four-body forces is considered to derive from the state dependence 
of the effective nucleon-nucleon interaction and the additional Pauli repulsion between more than two 
a-particles. However, they are short-range, and hence only act in compact configurations. The expec- 
tation values of those forces do not exceed 7 percent of the one of the corresponding two-body term, 
even for the ground state with the most compact structure, i.e. being the most sensitive to those forces. 

Employing the Gaussian expansion method [H] for the choice of variational basis functions, the 
total wave function of the 4a system is expanded in terms of Gaussian basis functions as follows: 

C,I/ 

<l>c{u) = S^[<^i,{ri,ui)^i^{r2,U2)]h2¥^is{r3,U3) ^, (23) 

where ri, r2 and are the Jacobi coordinates describing internal motions of the 4a system. S stands 
for the symmetrization operator acting on all a particles obeying Bose statistics, u denotes the set of 
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Table 2: The rms radii R and monopole transition matrix elements to the ground state M(EO) in 
units of fm and fm^, respectively. -Rexp. and M(EO)exp. are the corresponding experimental data. The 
finite-size effect of a particle is taken into account in R and M(EO) (see Ref. [23] for details). 
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size parameters ui, U2 and z/3 of the normalized Gaussian function, (pi{r, Ui) = Ni^i^.r'' exp {—Uir'^)Yim{r), 
and c the set of relative orbital angular momentum channels [[/i, /2]h25 ^3] J depending on either of the 
coordinate type oi K 01 H [H], where /i, I2 and Z3 are the orbital angular momenta with respect to the 
corresponding Jacobi coordinates. The coefficients A^^v) are determined according to the Rayleigh-Ritz 
variational principle. 

[MeV] 



Figure 12: Comparison of energy spectra between 
experiment and the present calculation. Two kinds 
of effective two-body nucleon-nucleon forces MHN 
and SW are adopted (see text). Dotted and dash- 
dotted lines denote the a±^^C and 4a thresholds, 
respectively. Experimental data are taken from 
Ref. |45], and from Ref. [46] for the 0^ state. The 
assignments with experiment are tentative, see, 
however, detailed discussion in the text. 



EXP. SW MHN 

Figure [T2] shows the energy spectrum with J'^ = O"*", which is obtained by diagonalizing the Hamil- 
tonian, Eq. (1211) . in a model space as large as given by 5120 Gaussian basis functions, Eq. (123|) (the 
other multipolarities, needing larger basis sets, are more difficult and shall be studied in future work). 
It is confirmed that all levels are well converged. With the above mentioned effective a-a forces, we can 
reproduce the full spectrum of 0^ states, and tentatively make a one-to-one correspondence of those 
states with the six lowest 0+ states of the experimental spectrum. In view of the complexity of the 
situation, the agreement is considered to be very satisfactory. We show in Table 121 the calculated rms 
radii and monopole matrix elements to the ground state, together with the corresponding experimental 
values. The M(EO) values for the 0^ and 0^ states are consistent with the corresponding experimental 
values. The consistency for the Oj]" state is within a factor of two. As mentioned above, the structures 
of the 0^ and Oj]" states are well established as having the a±^^C(0^) and a±^^C(2]^) cluster structures, 
respectively. These structures of the O2 and O^j" states are confirmed in the present calculation. We also 
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mention that the ground state is described as having a shell-model configuration within the present 
framework, the calculated rms value agreeing with the observed one (2.71 fm). 

On the contrary, the structures of the observed O4 , 0^ and 0^ states in Fig. [12] have, in the past, 
not clearly been understood, since they have never been discussed with the previous cluster model 
calculations [33 ESI EB] • Although Ref. [20], using the THSR wave function, predicts the 4a condensate 
state around the 4a threshold, it is not clear to which of those states it corresponds to. We will analyse 
the situation with the THSR wave function of [20] in a future publication |47j . 

As shown in Fig. [121 the present calculation succeeded, for the first time, to reproduce the O4 , 0^ 
and Og states, together with the 0^, 0^ and Ojj" states. This puts us in a favorable position to discuss 
the 4a condensate state, expected to exist around the 4a threshold. 

In Table 121 the largest rms value of about 5 fm is found for the Og state. Compared with the 
relatively smaller rms radii of the O4 and 0^ states, this large size suggests that the 0^ state may be 
composed of a weakly interacting gas of a particles [H] of the condensate type. 




20.0 



Figure 13: (Color online) The radial parts of 
single-a orbits with L = belonging to the largest 
occupation number, for the ground and Og states 
with MHN force. 




20.0 



Figure 14: (Color online) ry{r) defined by Eq. (l24l) 
for the Ot state with the MHN force. 



While a large size is generally necessary for forming an a condensate, the best way for its identi- 
fication is to investigate the single-a orbit and its occupation probability, which can be obtained by 
diagonalizing the one-body (a) density matrix as defined in [l9l [321 [231 EO]- As a result of the calcu- 
lation of the L = case, a large occupation probability of 61% of the lowest OS'-orbit is found for the 
0^ state, whereas the other five O"*" states all have appreciably smaller values, at most 25% (0^). The 
corresponding single-a 5* orbit is shown in Fig. [131 It has a strong spatially extended behaviour without 
any node (OS*). This indicates that a particles are condensed into the very dilute OS* single-a orbit, see 
also Ref. [51] . Thus, the 0^ state clearly has 4a condensate character. We should note that the orbit is 
very similar to the single-a orbit of the Hoyle state [321 [23] . We also show in Fig. [13] the single-a orbit 
for the ground state. It has maximum amplitude at around 3 fm and oscillations in the interior with 
two nodal (25*) behaviour, due to the Pauli principle and refiecting the shell-model configuration. 

In order to further analyze the obtained wave functions, we calculate an overlap amplitude, which 
is defined as follows: 

yir) = ( \^-^^^Y,if')^,C'C)] ^ ^(Og+)\. (24) 



J 



Here, $i(^^C) is the wave function of ^^C, given by the 3a OCM calculation [23], and r is the relative 
distance between the c.o.m. of ^^C and the a particle. From this quantity we can see how large is 
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the component in a certain a+^^C channel which is contained in our wave function fl23l) for Og'. The 
amphtudes for the 0^ state are shown in Fig. (TH It only has a large amplitude in the a+^^C(0^) 
channel, whereas the amplitudes in other channels are much suppressed. The amplitude in the Hoyle- 
state channel has no oscillations and a long tail stretches out to ~ 20 fm. This behaviour is very similar 
to that of the single-a orbit of the Ofj' state discussed above. 

The a decay width constitutes a very important information to identify the 0^ state from the 
experimental point of view. It can be estimated, based on the i?-matrix theory, with the overlap 
amplitude Eq. [S2]- We find that the total a decay width of the state is as small as 50 keV 
(experimental value: 166 keV). This means that the state can be observed as a quasi-stable state. 
Thus, the width, as well as the excitation energy, are consistent with the observed data. All the 
characteristics found from our OCM calculation, therefore, indicate that the calculated 6th O"*" state 
with 4 alpha condensate nature can probably be identified with the experimental 0^ state at 15.1 MeV. 

Finally we discuss the structures of the 04" and 0^ states. Our present calculations show that the Oj 
and 0^ states mainly have a+^'^G^Q^) structure with higher nodal behaviour and Q!-|-^^C(1~) structure, 
respectively. Further details will be given in forthcoming work. The calculated width of the 0^ is ~ 150 
keV, which is quite a bit larger than that found for the 0^ state ~ 50 keV. Both are qualitatively 
consistent with the corresponding experimental data, 600 keV and 185 keV, respectively. The reason 
why the width of the Oj state is larger than that of the 0^ state, though the 0^ state has lower excitation 
energy, is due to the fact that the former has a much larger component of the a+^^C(0]'') decay channel, 
reflecting the characteristic structure of the O4 state. The 4a condensate state, thus, should not be 
assigned to the O4 or 0^ state [53] but very likely to the 0^ state. 

In conclusion, the investigation of the 0"'"-spectrum with the 4a OCM calculation succeeded in 
describing the structure of the full observed 0"*" spectrum up to the 0^ state in ^^O. The 0"*" spectrum 
of ^^O up to about 15 MeV is now essentially understood, including the 4a condensate state. This is 
remarkable improvement concerning our knowledge of the structure of ^^O. We find that the 0^ state 
above the 4a threshold has a very large rms radius of about 5 fm and has a rather large occupation 
probability of 61% of four a particles sitting in a spatially extended single-a OS* orbit. The wave function 
has a large a+^^C amplitude only for ^^C*, i.e. the Hoyle state. These results are strong evidence of 
the Og state, which is a new theoretical prediction, for being the 4a condensate state, i.e. the analog 
to the Hoyle state in ^^C. Further experimental information is very much requested to confirm the 
novel interpretation of this state. Also independent theoretical calculations are strongly needed for 
confirmation of our results. 
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Figure 15: Alpha-particle mean-field potential for three a's in ^^C and six a's in ^^Mg. Note the lower 
Coulomb barrier for ^^Mg (from Ref. j56j). 
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In principle, one could go on, inreasing the number of a-particles: ^"Ne, ^^Mg, etc. However, one 
easily imagines that the complexity of the calculations quickly becomes prohibitive. In order to get 
a rough idea what happens for more a-particles, drastic approximations have to be performed. One 
such approximation is to consider the a-particles as ideal inert bosons and to treat them in mean field 
approximation. This then leads to the Gross-Pitaevskii Equation (GPE) which is widely employed 
in the physics of cold atoms [51] • One interesting question that can be asked in this connection is: 
How many a's can maximally exist in a self-bound a-gas state? Seeking an answer, we performed a 
schematic investigation using an effective a-a interaction within an a-gas mean-field calculation of the 
Gross-Pitaevskii type [55] • The parameters of the force were adjusted to reproduce our microscopic 
results for ^^C. The corresponding a mean-field potential is shown in Fig. [151 One sees the OS-state 
lying slightly above threshold but below the Coulomb barrier. As more a-particles are added, the 
Coulomb repulsion drives the loosely bound system of a-particles farther and farther apart, so that the 
Coulomb barrier fades away. According to our estimate [SB], a maximum of eight to ten a-particles 
can be held together in a condensate. However, there may be ways to lend additional stability to such 
systems. We know that in the case of *Be, adding one or two neutrons produces extra binding without 
seriously disturbing the pronounced a-cluster structure. Therefore, one has reason to speculate that 
adding a few of neutrons to a many-a state may stabilize the condensate. But again, state-of-the-art 
microscopic investigations are necessary before anything definite can be said about how extra neutrons 
will influence an a-particle condensate. A study in this direction is given in Ref. [ST] for ^^C. 



7 Alpha-particle condensates in finite nuclei and the alpha- 
particle occupation numbers 
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Figure 16: (Color online) Occupation of the 
single-a orbitals of the ground state of ^^C 
compared with the Hoyle state [23] . 




Figure 17: Momentum distribution of the a particle 
for the ground (solid line) and 0^ (dotted line) states. 



As already mentioned, constructing an a-particle density matrix p(-R, -R') by integrating out of 
the total density matrix all intrinsic a-particle coordinates and diagonalising the result, one finds that 
the corresponding 05" a-particle orbit in ^^C is occupied to more than 70 percent by the three a- 
particles [321 [23]. This is a huge percentage, affirming the almost ideal a-particle condensate nature of 
the Hoyle state. By contrast, even at zero temperature only ten percent of the particles in superfiuid 
^He belong to the condensate (which is nevertheless a macroscopic supply of condensed particles). To 
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add further perspective to the picture, in the ground state of ^^C the a-particle occupations are equally 
shared between S, D and G orbits, thus invalidating a condensate picture for the ground state. The 
occupation numbers for the ground and Hoyle states are shown in histogramm format in Fig. [121 The 
difference between the Hoyle state and the ground state is seen to be spectacular. In the Hoyle state 
the OS'-occupancy is at least an order of magnitude(!) higher than for any other orbit. This is one 
of the main typical features of Bose-Einstein condensation, even in strongly correlated Bose systems 
where there may be a strong depletion of the condensate, like in superfiuid ^He. On the other hand, 
the ground state occupancies can be explained quite well with the standard shell model [23]. It should 
also be noted that the ground state of ^^C is reasonably well reproduced by our theory (see Table [I]). 
A further strong indication of the condensate-like behavior of the a-particles in the Hoyle state is their 
momentum distribution, which is much narrower, almost delta-function-like, than in the ground state, 
see Fig. [T71 On the same line of investigation one finds, as already mentioned in section 6, that the 6-th 
0"*" state in ^^O has 61% OS occupancy of a-particles. 
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Figure 18: Radial parts of the sing le-a S orbits, (a) of the Hoyle state (^^c), (b) of the 0^ state in ^^O 
and of the ground states (c) in ^^C and (d) in ^^O. 



In Fig. [TSVa) and (b), we show radial parts of the single-a 5* orbits of the Hoyle state [23] and 
the 0^ state in ^^O [39], respectively. We see an almost identical shape! Of course, the extension is 
slightly different because of the smallness of the system. The nodeless character of the wave function is 
very pronounced and only some oscillations with small amplitude are present in ^^C, reflecting a weak 
influence of the Pauli principle between the a's! On the contrary, we show in Fig. [TST c) and (d), radial 
parts of the single-a S orbits of the ground states in ^^C [23] and ^^O [39], respectively. Due to its much 
reduced radius the "a-like" clusters strongly overlap, producing strong amplitude oscillations which 
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take care of antisymmetrisation between clusters. Again this example very impressively demonstrates 
the condensate nature of the Hoyle state and the 0^ state in ^^O. 

It is worth noting at this point that the definition of the a-particle density matrix, that is of a 
self bound system, is somewhat ambiguous and that different definitions may lead to different answers 
for the occupancies. This question has recently been debated in a number of papers [58l [59], [50] . 
Our conclusion concerning this point is that one should use for the definition of the internal density 
matrix an orthogonal systems of independent coordinates, like they are for instance given by the Jacobi 
coordinates. This definition insures the physically very reasonable boundary condition that, given a 
system exhibits an ideal Bose condensate in the laboratory system, it so remains in the internal system 
and both description become equivalent in the thermodynamic limit [50]. Jacobi coordinates have been 
used to evaluate the occupancies of a-particles in nuclei mentioned above [231 1221 [Ml [12] . 

8 Reduction of the a-condensate with increasing density 

The properties of a matter can be used to frame the discussion of the structure of na nuclei. As 
described in the preceding section, computational studies of these nuclei based on THSR cluster states 
have demonstrated that an a condensate is established at low nucleon density. More specifically, states 
lying near the threshold for decomposition into a particles, notably the ground state of ^Be, ^^C in 
the 0^ Hoyle state, and corresponding states in ^^O and other na nuclei are dilute^ being of low 
mean density and unusually extended for their mass numbers. We have shown quantitatively within 
a variational approach that a-like clusters are well formed, with the pair correlation function of a-like 
clusters predicting relatively large mean distances. For example, in determining the sizes of the ^^C 
nucleus in its 0^ (ground) state and in its 0^ excited state, we obtained rms radii of 2.44 fm and 3.83 
fm, respectively. The corresponding mean nucleon densities estimated from 36/47rrj?j^g are close to the 
nuclear-matter saturation density po = 0.16 nucleon/fm^ in the former state and 0.03 nucleon/fm^ in the 
latter. The expected low densities of putative alpha- condensate states are confirmed by experimental 
measurements of form factors |33j . 

All of our considerations indicate that quartetting is possible in the low-density regime of nucleonic 
matter, and that a condensates can survive until densities of about 0.03 nucleons/fm^ are reached. Here, 
we are in the region where the concept of a matter can reasonably be applied [GQl [61] . It is then clearly 
of interest to use this model to gain further insights into the formation of the condensate, and especially 
the reduction or suppression of the condensate due to repulsive interactions. We will show explicitly 
that in the model of a matter, as in our studies of finite nuclei, condensate formation is diminished 
with increasing density. Already within an a-matter model based on a simple a — a interaction, we can 
demonstrate that the condensate fraction - the fraction of particles in the condensate - is significantly 
reduced from unity at a density of 0.03 nucleon/fm'^ and essentially disappears approaching nuclear 
matter-saturation density. 

The quantum condensate formed by a homogeneous interacting boson system at zero temperature 
has been investigated in the classic 1956 paper of Penrose and Onsager [52] who characterize the 
phenomenon in terms of off-diagonal long-range order of the density matrix. Here we recall some 
of their results that are most relevant to our problem. Asymptotically, i.e., for \r — r'\ ~ oo, the 
nondiagonal density matrix in coordinate representation can be decomposed as 



In the limit, the second contribution on the right vanishes, and the first approaches the condensate 
fraction, formally defined by 




(25) 



Po 




(26) 
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Penrose and Onsager showed that in the case of a hard-core repulsion, the condensate fraction is 
determined by a filhng factor describing the ratio of the volume occupied by the hard spheres. They 
applied the theory to liquid ^He, and found that for a hard-sphere model of the atom-atom interaction 
yielding a filling factor of about 28%, the condensate fraction at zero temperature is reduced from unity 
(its value for the noninteracting system) to around 8%. (Remarkably, but to some extent fortuitously, 
this estimate is in rather good agreement with current experimental and theoretical values for the 
condensate fraction in liquid ^He.) 

To make a similar estimate of the condensate fraction for a matter, we follow Ref. [2] and assume an 
"excluded volume" for a particles of 20 fm^. At a nucleonic density of po/3, this corresponds to a filling 
factor of about 28%, the same as for liquid ^He. Thus, a substantial reduction of the condensate fraction 
from unity (for a noninteracting a-particle gas at zero temperature) is also expected in low-density a 
matter. 

Turning to a more systematic treatment, we proceed in much the same way as Clark and coworkers 
j60j . referring especially to the most recent study with M. T. Johnson. Adopting the a — a interaction 
potential 

K.(r) = 475 e-(°-^'^/^"^)'MeV - 130 e-(°-^^^"/^"^)'MeV (27) 

introduced by Ali and Bodmer [63], we calculate the reduction of the condensate fraction as function 
of density within what is now a rather standard variational approach. Alpha matter is described as an 
extended, uniform Bose system of interacting a particles, disregarding any change of the internal struc- 
ture of the a clusters with increasing density. In particular, the dissolution of bound states associated 
with Pauli blocking (Mott effect) is not taken into account in the present description. 

The simplest form of trial wave function incorporating the strong spatial correlations implied by the 
interaction potential (1271) is the familiar Jastrow choice. 



^{r,,...,rA) = l[fi\r.-r,\). (28) 



The normalization condition 



i<j 



oo 

47rpa / [fir)-l]r^dr = -l, (29) 







in which is the number density of a-particles, is imposed as a constraint on the variational wave 
function, in order to promote the convergence of the cluster expansion used to calculate the energy 
expectation value [M]. In the low-density limit, the energy functional [binding energy per a cluster as 
a functional of the correlation factor /(r)] is given by 



E[f] = 27rp, / 

^0 




df{ry' 



fir)V^ir)}r'dr, (30) 



dr 

where nia is the a-particle mass, while the condensate fraction is given by 

(31) 



po = exp |-47rpa [f{r) - 1]^ r^dr^ 



The variational two-body correlation factor / was taken as one of the forms employed by Clark and 
coworkers [60], namely 

/(r) = (1 - e^"'')(l + be-"' + ce"^"') . (32) 

At given density p, the expression for the energy expectation value is minimized with respect to the 
parameters a, b, and c, subject to the constraint (!29l) . It is important to note that these approximations, 
based on truncated cluster expansions, are reliable only at densities low enough that the length scale 
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associated with decay of /^(l) — 1 is sufficiently small compared to the average particle separation, 
which is inversely proportional to the cubic root of the density [BUI ED EH ES]- 

To give an example, for the nucleon density 4pa = 0.06 fm~^, a minimum of the energy expectation 
value ( |30l) was found at a = 0.616 fm~^, b = 1.221, and c = —5.306, with a corresponding energy per a 
cluster of —9.763 MeV and a condensate fraction of 0.750. The dependence of the condensate fraction 
on the nucleon density p = 4pa as determined in this exploratory calculation is displayed in Fig. 




nucleonic density [fm ] 



Figure 19: Reduction of condensate fraction 
in a matter with increasing nucleon den- 
sity. Exploratory calculations (full line) are 
compared with HNC calculations of Johnson 
and Clark [60] (crosses). For comparison, 
we show estimates of the condensate fraction 
in the 0^ (Hoyle) state of ^^C, according to 
Refs. [231 [32] (stars). 




Figure 20: Occupation of the 5* orbital as a function 
of density using the 3a OCM for [23]. 



The reduction of the condensate fraction of a matter to roughly 0.8 as given by our calculation at 
nucleonic density 0.03 fm~^ agrees well with results of Suzuki [32] and Yamada [23] for ^^C in the Hoyle 
0^ state. Using many-particle approaches to the ground-state wave function and to the THSR (0^) 
state of ^^C, the occupation of the inferred natural a orbitals is found to be quite different in the two 
cases. Roughly 1/3 shares (approaching equipartition) are found for the S, D, and, G orbits in the 
ground (0^*^) state, with a-cluster occupations of 1.07, 1.07, and 0.82, respectively. On the other hand, 
in the Hoyle (Oj) state, one sees enhanced occupation (2.38) of the S orbit and reduced occupation 
(0.29, 0.16, respectively) of the D and G orbits. This corresponds to an enhancement of about 70% 
compared with equipartition. 

To get a more extended analysis, OCM calculations have been performed [23] for studying the density 
dependence of the S'-orbit occupancy in the Hoyle state on the different densities p/po ~ {R{Of)^^p/R)^, 
in which the rms radius (R) of ^^C is taken as a parameter and R{0f)^^p=2.56 fm. A Pauli-principle 
respected OCM basis *°+™(z/) with a size parameter u is used, in which the value of u is chosen to 
reproduce a given rms radius R of ^^C, and the a density matrix p{r,r') with respect to \l/^-^^(z/) is 
diagonalized to obtain the 5'-orbit occupancy in the O"*" wave function. The results are shown in Fig. [201 
The S'-orbit occupancy is 70 ~ 80 % around p/po ~ (-R(Oi")exp/-R(0^)THSR)'^ ~ 0.21, while it decreases 
with increasing p/po and amounts to about 30 ~ 40 % in the saturation density region. Figure EH 
shows the radial behaviours of the S'-orbit with given densities. A smooth transition of the S-orbit 
is observed, with decreasing p/po, from a two-node S'-wave nature (p/po ~ 1.18) in Fig. [2Tl(a) to the 
zero-node S'-wave one {p/po — 0.15) in Fig. [2]T d) [23] . The feature of the decrease of the enhanced 
occupation of the S orbit is in striking correspondence with the density dependence of the condensate 
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fraction calculated for nuclear matter (see Fig. [T9|) . 




Figure 21: Radial behaviors of the S orbit in the ^^0(0+) state with (a) R = 2.42 fm (p/po ~ 1-18), (b) 
R = 2.70 fm (p/po ~ 0.85), (c) R = 3.11 fm (p/po ~ 0.56), and (d) R = 4.84 fm (p/po ~ 0.15), where 
R denotes the nuclear radius of the ^^C(0"^) state. 



A more accurate and reliable variational description of a matter can be realized within the hypernetted- 
chain (HNC) approach to evaluation of correlated integrals; this approach [601 largely overcomes 
the limitations of the cluster-expansion treatment, including the need for an explicit normalization 
constraint. Such an improved approach is certainly required near the saturation density of nuclear 
matter, where it predicts only a small condensate fraction [60]. Of course, at high densities the simple 
Ali-Bodmer interaction [63] ceases to be valid, and it becomes crucial to include the effects of Pauli 
blocking. Once again, this conclusion reinforces the point that we can expect signatures of an a con- 
densate only for dilute nuclei near the threshold of na decay, but no signatures from configurations 
with saturated density. 



9 Conclusions, Discussion, Outlook 

We have investigated the role that pairing and multiparticle correlations may play in nuclear matter 
existing in dense astrophysical objects and in finite nuclei. A complete and quantitative description 
of nuclear matter must allow for the presence of clusters of nucleons, bound or metastable, possibly 
forming a quantum condensate. In particular, quartetting correlations, responsible for the emergence of 
a-like clusters, are identified as uniquely important in determining the behavior of nuclear matter in the 
limiting regime of low density and low temperature. We have calculated the transition temperature for 
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the onset of quantum condensates made up of a-like and deuteron-like bosonic clusters, and considered 
in considerable detail the intriguing example of Bose-Einstein condensation of a particles. It turns out 
that contrary to pairing, quartet condensation primarily exists in the BEC phase at low density. In 
which way quartet condensation is lost by increasing the density is still an open question. However, it is 
clear that there can not exist a condensate of quartets with a long coherence length for arbitrarily small 
attraction as this is the case for pairing in the BCS phase. It is inevitable that under increasing density 
or pressure, the bound a, d, or other nuclidic clusters present at low density experience significant 
modification due to the background medium (and eventually merge with it). We have shown how self- 
energy corrections and Pauli blocking alter the properties of cluster states, and we have formulated a 
cluster mean-field approximation to provide an initial description of this process. One result of special 
interest is the suppression of the a-like condensate, which is dominant at lower densities, as the density 
reaches and exceeds the Mott value, allowing the pairing transition to occur. Even at lower densities 
a-particle condensation may be influenced by neutron excess, i.e. in the case of asymmetric nuclear 
matter. The study of a-particle condensation as a function of asymmetry remains a task for the future. 

A truly remarkable manifestation of a-particle condensation seems to be present in finite nuclei. 
Indeed, the so-called Hoyle state (0^) in at 7.654 MeV is very likely a dilute gas of three a-particles, 
held together only by the Coulomb barrier. This view is encouraged by the fact that we can explain 
all the experimental data in terms of a conceptually simple wave function of the quartet-condensate 
type. Within the same model, we also systematically predict such states in heavier n a nuclei, and the 
search is on for their experimental identification. In a recent study with OCM (Orthogonality Condition 
Model) we predicted the 6-th O"*" state of ^^O to be a strong candidate for a loosely bound four a-particle 
state [39]. The results of that study are presented and discussed here. The condensate feature of the 
Hoyle state in ^^C and Hoyle- like state in ^^O is born out by the calculation of the bosonic occupation 
numbers in diagonalising the bosonic density matrix. It is shown that the occupation of the OS* state 
of the a-particles is over 70% for the Hoyle state in ^^C and over 60% for the Hoyle-like state in ^^O. 

It is quite natural that such loosely bound a-particle states should exist up to some maximum 
number of a particles. We estimate that the phenomenon will terminate at about eight to ten a's as 
the confining Coulomb barrier fades away. However, there is the possibility that larger condensates 
could be stabilized by addition of a few neutrons. Indeed, consider ^Be, which, contrary to ^Be, is 
bound by ~ 1.5 MeV, still showing a pronounced two a-structure similar to the one of Fig. |4] (b). One 
could imagine ten a's or more, stabilised by two or four extra neutrons in a low density phase. However, 
even without being stabilised, if a compressed hot nuclear blob as e.g. produced in a central Heavy 
Ion collision expands and cools, it may turn on its way out, at a certain low density, into an expanding 
a condensed state where all a's are in relative S'-waves. This would be an analogous situation to an 
expanding Bose condensate of atoms, once the trapping potential has been switched off. Forthcoming 
analysis of dedicated experiments with high resolution multiparticle detectors like CHIMERA at LNS, 
Catania, will tell whether such scenarios can be realised [66]. Other possibilities of loose a-gas states 
may exist on top of particularly stable cores, like ^^0 or ^°Ca. Indeed in adding a's to e.g. ^"Ca, one will 
reach the a-particle drip line. Compound states of heavy N = Z nuclei of this kind may be produced in 
heavy ion reactions and an enhanced a-decay rate may reveal the existence of an a-particle condensate. 
Ideas of this type are presently promoted by von Oertzen [671 EH], and also M. Brenner [69], and A. 
Ogloblin [To]. However, coincidence measurements of multiple a's of decaying lighter nuclei like ^^O 
may also be very useful [7T] [72] to detect at least one additional a-condensate state beyond the only 
one that has been identified so far, namely the 0^-state in ^^C. 

In finite nuclei a-particle condensation has to be understood in the same sense as we say nuclei are 
superfluid inspite of the presence of only a limited number of Cooper pairs. However, in compact stars, 
formed by Supernova Explosions, macroscopic condensates of a-particles may be present. In the study 
of reference [2], a-particle phases are predicted for temperatures which can easily be below the critical 
temperatures obtained by our calculation displayed in Figs. [2] and [31 For the study of such macroscopic 
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condensates, the use of the Gross-Pitaevskii equation for interacting ideal bosons can be useful [56] but 
may be hampered by our poor knowledge of the density dependence of the effective a-a interaction. 
A genuine microscopic approach for a-particle condensation in infinite matter is demanded, since our 
wave function, Eqs. ([8]) and ([9]), only can handle a limited number of a's due to the necessary explicit 
antisymmetrisation. A preliminary route to this aim is outlined in this presentation. 

Another issue which may be raised in the context of a-particle condensation is the question, also 
discussed in condensed matter physics [73], whether a's condense as singles or as doubles, i.e. as ^Be? In 
microscopic studies of ^^C one, indeed, can see that in the Oj-state two of the three a's are slightly more 
closer to one another than to the third one [7^ [75l [76] . Without Coulomb repulsion ^Be may be bound 
by 2-3 MeV, however, the Coulomb repulsion makes it very slightly unbound. The question is definitely 
very interesting and deserves future investigation. However, quantitatively, if at all, a-a correlations 
constitute certainly only a slight modification over the present formulation of a-condensation. One may 
even ask whether these ^Be-like correlations are not an artifact of the calculation because in the RGM 
calculation of [15] ®Be-like correlations are definitely allowed, however, as already mentioned, our pure 
a-particle condensate wave function turned out to have almost 100 percent overlap with the one of the 
RGM calculation. 

As already outlined earlier, nuclear systems exhibit especially strong cluster and few body effects. 
This stems primarily from the fact that there are four different fermions with in addition more or less 
equal attraction among one another. However, this situation is not necessarily unique. In the past there 
have been speculations that in semi-conductors bi-excitons may condense rather than excitons [77]. A 
very promising field in this respect may be the possibility to trap in magneto-optical devices fermionic 
atoms in four different magnetic sub-states [3]. If in addition the interaction between all four fermions 
could be triggered, eventually with the help of Feshbach resonances, so that they are all attractive, then 
investigation of quartetting could become a domain of research as rich as is presently the investigation 
of pairing. Clearly, the study of quartet condensation only is at its beginning and perspectives are rich 
and manyfold. 

In conclusion, we see that the idea of a-particle condensation in nuclei has already triggered many 
new ideas and calculations, in spite of the fact that, so far, a compelling case for such a state has only 
been made in ^^C. Even so, the possible existence of a completely new nuclear phase in which a-particles 
play the role of quasi-elementary constituents is surely fascinating. Hopefully, many more a-particle 
states of nuclei will be detected in the near future, bringing deeper insights into the role of clustering 
and quantum condensates in systems of strongly interacting fermions. 
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